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Abstract. - We argue that Poisson statistics in logarithmic time provides an idealized descrip- 
tion of non-equilibrium configurational rearrangements in aging glassy systems. The description 
puts stringent requirements on the geometry of the metastable attractors visited at age t^. An- 
alytical implications for the residence time distributions as function of and the correlation 
functions are derived. These are verified by extensive numerical studies of short range Ising 
spin glasses. 



Introduction. - Aging is a key aspect of the dynamics in a host of complex glassy sys- 
tems IPI^llSlEEliniEIIHI ■ Its main characteristics are a slow, sub- linear drift of macroscopic 
averages with the age t^ of the system, i.e. the time elapsed after an initial thermal quench, 
and a concomitant age dependence of the correlation and response functions. The latter shows 
that a quasi-stationary fluctuation regime for short observation times t < <„, is followed by 
off-equilibrium dynamics for t > ty^. This behavior can be attributed to an increasing degree 
of thermal stability of the dynamically relevant attractors or landscape valleys as the age tw 

grows inmniiniini. 

Within each valley, the pseudo-equilibrium fluctuations are determined by the free energy. 
However, the inter- valley non-equilibrium dynamics could be practically irreversible, and hence 
qualitatively different. Adopting the simplified view that the selection is truly irreversible 
allows us to formulate a theoretical description of large non-equilibrium rearrangements, or 
'quakes', which borrows from studies of non-thermal models 13 14 and which supplies a basis 
to the striking similarity of complex dynamics seen in very different systems with multiple 
metastable states. The key feature is dynamical invariance to shifts in the logarithm of time. 
This implies that the time spent in thermally metastable attractors is power-law distributed 
with an exponent a. Furthermore, the distribution has a t/tw aging behavior, often called full 
or pure aging. The exponent a depends on the system size but not on the temperature. The 
further assumption that the memory of the initial configuration decays as exp{—^k) with the 
number k of quakes leads to a power-law decay and a,t/tw form of the correlation function 
for t > t^. The corresponding non-equilibrium exponent A is jointly determined by a and 7. 

These theoretical results are checked through extensive numerical simulations of short 
ranged spin glass systems. While direct experimental evidence is not yet available, possibly 
related large-scale irreversible rearrangements similar to quakes have recently been identified 
in aging colloidal gels (3- 
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Basic formalism. - Beyond its original macro-evolutionary context, the concept of 'punc- 
tuated equilibrium' |15| is applicable to any stochastic dynamics where long periods of stag- 
nation are dotted by fast irreversible changes. This occurs e.g. in micro-evolution "16" 17' and 
in certain driven dissipative models |18II14| . Within a population evolving in a fitness land- 
scape |17j , fitness records achieved by an individual through random mutations are amplified 
to the macroscopic level by darwinistic selection and trigger punctuations. Since the same 
statistics describes both punctuations and records, the analytical description available for the 
latter |13| can be applied to the former. Below, we discuss how thermal noise may similarly 
induce large and irreversible quakes in glassy systems. 

To summarize the basic ideas of record statistics IS , consider times > I and i > 0, 
and independently draw for each integer time a random number from any fixed distribution 
not supported on a finite set. The attempts whose outcome exceeds all previous outcomes are 
marked as records. For large times, the probability that k such events occur in [i^,, i^, -I- t] is 
then 



where a — I for a single sequence of attempts. The total number of records produced by n in- 
dependent sequences is a sum of log-Poisson variables, and hence itself log-Poisson distributed, 
with a = n. 

Equation can be recognized as a Poisson distribution with log{t^+t) — \og{t^) replacing 
the linear difference t between the time arguments, or, for short, a log-Poisson distribution(^). 
Equation implies that the rate of events decays as a/(t -f t^,), while the logarithmic rate 
of events approaches the constant value a for t 3> . 

Glossing over the discreteness of the time arguments, we can now pick thermal noise as the 
(arbitrary) source of randomness producing the records. Since landscape induced correlations 
are negligible if records are intercalated by long intervals spent dwelling near a local energy 
minimum, the number of energy records observed in a fixed time interval is expected to be 
proportional to a log-Poisson variable ^Hl • The same applies if energies are measured relative 
to the current lowest energy minimum, producing the energy difference henceforth referred 
to as a 'barrier'. With this definition, not all barrier crossings carry physical significance. 
However, if quakes are log-Poisson distributed, they are naturally associated with the crossing 
of barriers of record magnitude. Two special landscape features sketched in Fig^would allow 
barrier records to induce quakes: i) The stability of the attractors successively selected, as 
gauged by e.g. an exit barrier, must only increase marginally |18II14| . i.e. in principle by an 
infinitesimal amount. This requires that shallow attractors be predominant, and implies that 
the actual magnitude of the records generating the quakes is immaterial(^). ii) A large and 
negative energy change associated to a quake, ensures its irreversibility. Such changes are 
observed in both atomistic glass models [201 and spin glasses |18j . 

The only adjustable quantities in this landscape cartoon are the parameter a and the form 
of the quake induced decorrelation. The noise independence of the statistics restricts a to 
be a temperature independent quantity. An interesting increase with system size can arise in 
spatially extended systems, where sufficiently distant regions can evolve independently. 

In the sequel, we first deal with some theoretical consequences of the marginal stability of 
the attractors and irreversibility of the quakes. To make contact with experiments, we then 

{l)The equat ion was derived in 1131 for tw = 1, but the generalization to the propagator Pn{tw,tw + t) is 
straightforward and completely analogous to the standard case. 

(^)The predominance of shallow attractors concurs with the observation that thermal quenches produce 
shallow minima in glassy systems, and fits the anomalous sensitivity of complex relaxation to small mechanical 
or thermal perturbations implied by the so-called reset or memory effects 1191 . 
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Fig. 1 Fig. 2 

Fig. 1 - An idealized landscape where the attractors are symbolized by wedges. The arrows represent 
'quakes' leading from one attractor to the next. They entail a large energy change — the vertical 
displacement — and are therefore irreversible. The log-Poisson statistics of the quakes originates from 
the infinitesimal increase of the metastability of the attractors selected, which is rendered by a slight 
increase in the bottom to top distance (depth) of each trap. 

Fig. 2 - The shown curves all describe the distribution of a suitably rescaled residence time in the 
attractors — or valleys — of the sort sketched in Fig. The circles are the prediction of the log-Poisson 
theory — Eqs. Q and ^ — and correspond to full aging (/x = 1). The empirical data for the spin 
glass landscape are plotted as lines, and describe the scaled residence time tr/^m in valleys entered at 
10^ < tw < 10^ . They merge into a seemingly continuous line due to the quality of the data collapse. 
The insert shows the same data with the less satisfactory p = 1 scaling. 

summarize an empirical way to identify quakes in glassy systems |18| and apply it to spin 
glasses. 

Attractor residence times. - The distribution of the residence or trapping time tr in 
valleys of an energy landscape has been considered theoretically |11I12| . but has not previously 
been measured. 

Let tk mark the occurrence of the fc'th event in a log-Poisson process with logarithmic rate 
a. In perfect analogy to the usual Poisson process, the 'log-waiting time' log(tfc) — log(ife_i) 
has the exponential distribution 

Prob(log(tfe/tfc_i) < x) = 1 ~exp{-ax), 

which is mathematically equivalent to the event distribution given in Eq. ^ . We now consider 
the distribution of the residence time tr — tk — ifc-i for traps entered at tk-i — t-^j. As 
tr/tw = tk/tk-1 - 1, we find 

Proh{tr/t^ <x)^l-{x + l)-", (2) 

which corresponds to full aging behavior. Equation. with a = 2.3 is plotted (circles) in 
the main panel of Fig. [3 together with spin glass data (lines). These are scaled in a slightly 
different manner, to account for a small but systematic deviation from full aging. 

The power-law form of Eq. ||2J) was suggested by Bouchaud , with the added and crucial 
assumption of an infinite average residence time. In our case, the average residence time (tr) 
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is equal to tu,/(a — 1) which is finite and w for the observed a « 2. As a consistency check, 
we note that the largest relaxation time r^q in a valley must be smaller than {tr)^ for local 
equilibration to apply. When the latter is of order i^,, Teq < follows. In the lack of special 
symmetries |21| . Teq = Tcorr, where Tcorr characterizes the short time decay of correlation 
functions, either as a cut-off parameter or directly as the time constant of an exponential 
decay. Hence Tcorr < tw as expected. 

Scaling of correlations. - Additional information on the effect of the quakes enters the 
calculation of thermally averaged correlations functions. Let c(m; m+k) be the configurational 
overlap between the lowest energy configurations in valley m and m + k. With probability 
Pm{^,tyj)Pk{tw,tw + t), m and m + k are the 'current' attractors at times and + t, 
whence the non-equilibrium part of the configurational autocorrelation function is given by 

C(t^,t^ +t)=^ ^c{m;m + k)P,n{l,t^)Pk{t^,t^ +t). (3) 

m — 1 k—Q 

Of special interest is c(l; 1 -I- k), the average overlap between the first and the (1 -I- fc)'th 
valley. This is the sole contribution to Eq. Q for t^, = 1, in which case C(l, 1 + t) additionally 
describes the magnetization decay of the fully polarized configuration of a spin glass ■ The 
exponential form 

c(l; l + k)^{l- Coo)e-T(^)'= + Coo (4) 

anticipates the results in Fig. ^ An exponential decay means that each event induces the 
same relative change, and leads to 

oo 

1 + t) - ^ c(l; 1 + fc)Pfe(l, l + t) = {l- Coo)t-^ + Coo, (5) 

k=Q 

where the non-equilibrium exponent A(T) is given by 

A(T) =a(l-e-^(^)) «a7(T). (6) 

Equation @ thus links a dynamical exponent to the landscape geometry. The second approx- 
imate equality given in Eq. JHl is only justified for 7 <^ 1. 

The further assumption that c{m; m + k) is altogether independent of m amounts to 
translational invariance with respect to m, and hence with respect to additive shifts in the 
logarithm of ty^. These are equivalent to a multiplicative rescaling of t^. Not surprisingly, 
Eq. Q yields the scale invariant and full aging form 

C{t^,t^ + t)o,(^^±^^ \ (7) 

Landscape exploration. - How should valleys and quakes be concretely defined given a 
time sequence of data? One approach is the protocol developed [18, for numerical data describ- 
ing isothermal relaxation after a deep quench. The basic observation is that at sufficiently low 
T pseudo-equilibrium fluctuations repeatedly visit the same local energy minimum. Hence, 
the attainment of an energy value lower than all previous values implies a non-equilibrium 
event, i.e. an ongoing quake. Several closely spaced minima records will often occur as a quake 
slides downhill among which the one preceding a barrier record is likely to be the most phys- 
ically significant, as it — ideally — coincides with a new local minimum. By the same token, in 
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a series of barrier records corresponding to an uphill climb, the last record preceding a low 
energy record is most significant, as it marks the access to a new local minimum, and hence 
the start of a quake. The times at which the above types of event occur provide an empirical 
definition of both valleys and quakes. 

This procedure yields a null result if i) only a single energy minimum is present; or ii) the 
initial condition is the ground state irrespective of the number of local minima present; or iii) 
all the minima and/or barriers are equal. Also, iv) if the trajectories can bypass the energy 
barriers and edge their way into states of lower energy through 'entropic barriers', the number 
of valleys discovered remains less or equal to one at all times. 

The actual outcome is very different, as illustrated in Fig.|21 where the estimated average 
and variance of the number nv of valleys visited is plotted versus the logarithm of time. An 
equality between variance and average is a key feature of Poisson distributions. If we neglect 
the initial curvature, the estimated quantities are very nearly proportional to logt, indicating 
that nv itself is proportional to a log-Poisson variable through a constant of order one. We 
emphasize that log-Poisson statistics is an idealized low T description, and that systematic 
corrections in the form of a decaying logarithmic rate of events must be expected, since each 
valley has a growing probability to contain the true ground state. Since a is not truly constant, 
the single value required by the log-Poisson theory cannot be estimated accurately. 

As T grows, local thermal equilibration within each attractor is not granted, and low en- 
ergy records lose connection to landscape minima. A check of the energy and configurational 
changes associated with the quakes, of the sort performed in |18| is required to indepen- 
dently establish their physical relevance. We additionally checked that the low energy record 
configurations are close (albeit usually not identical) to actual valley minima. 

Numerical results. - We performed extensive simulations of short-ranged Ising spin 
glasses on cubic lattices {L^ spins) with periodic boundary conditions and symmetric Gaus- 
sian couplings with unit variance, see also ^S]. The Waiting Time Method with single 
spin flip dynamics, which is dynamically equivalent to the Metropolis algorithm but much 
faster at low T, was used. The 'intrinsic' and size independent time variable of this method 
corresponds to the number of Monte Carlo (lattice) sweeps in the Metropolis as well as to the 
physical time of a real experiment. In all our runs, we skip the data within the first 10 time 
units in order to let the system settle down from its random T = oo initial configuration. 

The set of data shown in Fig. O was obtained for a system of linear size L = 30 at 
temperature T = 0.5. Numerous similar plots (not shown) were obtained for temperatures in 
the range 0.2 < T < 0.8, for 8 < L < 30 and < 10*^. In most cases 20000 independent 
runs were performed for each value of T and L, removing any visible statistical flutter. When 
a is calculated as the logarithmic slope of the number of valleys seen on average in the time 
window [10*, 10^], its size dependence for L — 13, 16,20 and 30 is found to be a = 1.17(1), 
1.37(1), 1.58(1) and 1.63(1), respectively. While the growth clearly tapers off, we cannot 
with confidence conclude on the true asymptotic behavior for large L. We note, however, 
that A(T) « aj{T) should remain finite in the thermodynamic limit. Hence, the growth of 
a must either remain bounded, or be compensated by a corresponding decrease of "f(T). To 
see that the latter is likely the case, we note that 7 — — log(l — 2Hii/L'^), where is the 
Hamming distance between the local minima at the bottom of consecutive {k = 1) valleys. 
In 3 d, Hd oc L^ ® was found in TS*, whence, using Hd/L'^ << 1 one obtains the slow decay 
7OCL-0-2. 

In order to obtain a proper log-Poisson variable, the number of valleys must be divided by 
~ 0.7, the ratio between variance and average, and the empirical value of a must be similarly 
rescaled to the effective value a ~ 2.0 for a 16"^ system, and a « 2.3 for a 30"^ system. 
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Fig. 3 - The average and variance of the number of vaUeys nv visited in time [10,tuj]. Their ratio 
is shown as a function of time in the insert. After a short transient, average and variance are nearly 
proportional to one another and to log , indicating that ny is proportional to a log-Poisson variable. 
In the time interval [10*, 10®], the logarithmic slope is q = 1.63(1). 

Fig. 4 - The main panel shows that the assumed exponential form of (c(fc) — c(cx>))/ (1 — c(c>o)) versus 
k is supported by the data, for temperatures T = 0.30 (upper curve) 0.35, 0.40 . . . 0.70. The lines are 
obtained by a two-parameter fits, and the symbols correspond to data points. The insert shows the 
temperature dependence of the non-equilibrium exponent A(r, = 10) according to 22 (circles) 
and as given by Eq. © with a — 2.0. The age argument — 10 was chosen for the comparison, as 
it corresponds to the earliest time at which our data were collected. All temperatures are divided by 
To = 0.9. 



Next, we analyze the distribution of the residence time tr within valleys entered at time 
[0.95tu,, l.OSiu,] for 10^ < < 10^. The observed frequencies of the scaled quantities tr/tl^ 
with /i = 1.05 are obtained from the same runs as in Fig. 13 They are plotted as lines in 
the main panel of Fig. |21 and the excellent quality of the resulting collapse gives the false 
impression of a single continuous line. The same distributions are also plotted in the insert, 
this time with ^ = 1, which is the theoretical prediction of the log-Poisson statistics. We see 
that while fJ, = 1 fails to produce a complete data collapse, it nevertheless accounts for the 
main features of the distribution. The corresponding theoretical prediction of Eq. ^ with 
a = 2.3 is given by the circles in the main panel. 

We finally consider the correlation decay. The overlap c(l; fc-l- 1) — for short c(fc) — between 
the lowest lying configurations of the first and the (1-1- fc)'th valley was calculated separately 
for temperatures T = 0.3, 0.35, . . . , 0.7 starting in each case from a random initial condition. 
In this case, data are for each T averages over 2000 trajectories belonging to 20 different 
realizations of the couplings. Since no data were taken for the first 10 time steps, the smallest 
age value considered is = 10, rather than the theoretical value t^, = 1. A two parameter 
non-linear fit to the exponential form in Eq. Q was performed, producing the lines shown in 
the main panel of Fig. ^ while the actual data are given by the symbols. The main error in 
fitting c(fc) is systematic and originates from the poor sampling of late events, which is dealt 
with by introducing a cut-off at fc = 12. The fitted values of 7 and Coo vary a few percent as 
the position of the cut-off is varied slightly around this value. The insert of Fig. 0] compares 
the temperature dependence of the predicted non-equilibrium exponent A w aiiT) (squares) 
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with the corresponding result by Kisker et al. j22| . The comparison is primarily qualitative, 
since we are unable to determine a with great precision. We used a = 2.0 in this plot. 

The (implied) temperature dependence of 7(T) shows that, for the same number of quakes, 
configurations decorrelate faster at higher temperatures, as one would expect. We note in 
passing that Coo systematically decreases with T, and that Coo would equal the Edwards- 
Anderson order parameter if the simple exponential decay of c(fc) were to continue indefinitely. 

While the exponential decorrelation form given by Eq. |(7J) — and the self-similar picture of 
the landscape it implies — is required for full aging, it is is not likely to be exactly fulfilled on 
all time scales. Yet it agrees with the full aging behavior observed in spin glasses which 
highlights the dynamical nature of their behavior. Interestingly, strong superaging behavior is 
found in gels Hence, either the latter systems are closer to thermal equilibrium than spin 
glasses, or considerable deviations from the marginal stability increase of the visited attractors 
occur. 

* * * 
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